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TECHNICAL. NOTE 3897 


INCOMPLETE TIME RESPONSE TO A UNIT IMPULSE AND ITS 
APPLICATION TO LIGHTLY DAMPED LINEAR SYSTEMS 
By Janies J. Donegan and Carl R. Hubs 


SUMMARY 


A study Is made of the use of the incomplete time response to a unit 
impulse as a device- to analyze lightly damped linear systems . It was 
found that the incomplete time response to a unit impulse computed from 
an incomplete time response and complete input adequately describes a 
lightly damped linear system without requiring an estimated end correc- 
tion. The incomplete response to a unit impulse may be used with 
Duhamel's integral to predict the response of the system to an arbitrary 
input up to at least 90 percent of the cutoff time of the incomplete 
response. 


INTRODUCTION 


One result frequently desired from the analysis of flight time- 
history data is the frequency response of the aircraft. Some of the 
methods available for accurately obtaining this result are given in ref- 
erences 1 and 2. As shown In references 2 and 3, the information con- 
tained In the frequency response can be converted to the time plane in 
the form of the complete time response to a unit impulse without requiring 
previous knowledge of the transfer function relating the Input and out- 
put. This information Is valuable in either form (time or frequency) 
Inasmuch as it completely defines the system and thereby permits pre- 
dicting the system response to any arbitrary input. 

In some cases, however, an accurate frequency response is difficult 
to obtain. Such a case is often experienced when aircraft motions are 
so lightly damped that it is impractical to continue the recording pro- 
cess. until the motion has subsided. The analysis of such systems by 
Fourier methods requires an estimated end correction to the incomplete- 
time-response data in order to obtain the system frequency response. 
Inaccuracies in the estimation of this end correction may introduce signi- 
ficant errors in the computed frequency response and, hence, in the com- 
puted response to an arbitrary input. 
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For a lightly damped higher order linear system where it is impos- 
sible to obtain an accurate frequency response because of limited time 
data, methods of analysis are needed to define the system by using input 
and output data, It“wcruld be desired that such a method use the data as 
recorded and bypass the need for estimating end corrections. 

The primary purpose of this report is to describe one such analysis 
technique which is called the incomplete time response- to a unit impulse, 
"incomplete" signifying that the response to the- unit impulse is obtained 
from an incomplete time response. 
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SYMBOLS 

Fourier transform of h(t) 
time response to unit-impulse 

incomplete time response to. unit- impulse cut off at 
A seconds 

time response to unit impulse in sideslip 

real component of Fourier transform 
time, sec 
output quantity 

Fourier transform, / x(t)e” :I - a; *dt 

response of system to arbitrary input 
sideslip angle, deg 

sideslip angle, as function of time - 
input quantity 

Fourier transform of input 
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S(t) 


time history of input to linear system 

s r 


rudder angle, deg 

6r(t) 


time history of rudder input 

A 


cutoff time of time history, sec 

0 


phase angle, deg 

0x5 


phase angle between output x and input 5 of frequency 
response; negative phase angles indicate lag 

0p6 r 


phase angle between output p and input 8 r 

CJU 


circular frequency, radians/sec 

1 x 1 
1 & r 

1-2-1 

M 

amplitude ratios 


A dot over a symbol denotes the derivative with respect to time. 
The absolute value of a term is denoted by | J . 

METHOD AND RESULTS 


In order to demonstrate the range of usef ul n ess of the incomplete 
time response to a unit impulse and to indicate the importance of end 
corrections , Fourier methods are first applied to a lightly damped mathe- 
matical system to show the effect on the results of cutting off the time 
histories at various points . The computed ratios of the Fourier trans- 
forms of output and input are compared as well as the time responses to 
a unit impulse. The same methods are then applied to a lightly damped 
response in sideslip as measured in flight, and the incomplete time 
response to a unit impulse is obtained. This incomplete time response 
is then used with the actual input to compute a time response which is 
compared with the measured aircraft time response. 

The Fourier transform, as used in this paper, is given as follows: 

x(ioj) = F x(t)e -lmfc dt + r x(t)e“ ic0fc dt 
Jq J A 
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where the integral with the limits' 0 to 7s Is^fche incomplete Fourier 
transform and the integral with the limits 7s to 00 is the Fourier 
transform of the end correction. 


Lightly Damped Mathematical System 

Figure 1 shows the exact calculated response of a lightly damped 
second-order system to a triangular input. This mathematical system is 
defined hy the following equation: 


x(t) + 0.2x(t) + 1.86x(t) = 0.75&(t) 


with the initial conditions x(0) =0 and x(o) =0. The exact frequency 
response obtained analytically for this system is given in figure 2 by 
the solid line which corresponds to a cutoff time An = <». This response 

hereinafter is called case 1. Two other computations representing the 
Fourier transforms x(iml 

are carried out by using the method of ref- 

5 ( io) 

erence 2 which is applied to a complete time history of the input and 
an incomplete time history of the output - of figure 1. In one case, the 
output was cut off at A 2 = 13-9 seconds and in the other, at 

7\j = 4.72 seconds'; these cases henceforth are referred to as case 2 
and case 3> respectively. In figure 2(a) the amplitude ratio resulting 
from the ratio of the transforms ^ ^ j is shown plotted against fre- 
quency, and in figure 2(b) the phase angle resulting from the ratio of 
the transforms is shown. The difference between cases 2 and 3 (the 
response cut off at . 7^ = 13.9 seconds and Aj =-4.72 seconds, respec- 
tively) and case 1 (the complete system A^ = °°) indicates the contri- 
bution of the Fourier transform of- the end corrections of cases 2 and 3 
to the complete -system frequency response. 


In order to compute the time response to a unit impulse by the 

method of reference 2, the real part of the Fourier transform Rj^-(io))J 

must be computed by multiplying the amplitude ratio shown in figure- 2(a) 
by the cosine of- the phase angle given in figure 2(b). The result of - 
this computation is shown in figure 3 for all three cases. A staircase 
function using a frequency interval Aco = 0.10 radian/sec was fitted 

to R|j~.(io)j for cases 2 and 3, and table I of reference 2 was used to 
compute the response to a unit impulse. 
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The resulting incomplete time responses to a unit impulse are shown 
in figure 4 along with the exact time response for the complete system. 
Within the cutoff interval the computations could have "been made more 
exact by taking a smaller interval Am in fitting the staircase function 

to R|j^(iasjJ . The largest error in h(t) occurs at t = 0 because 

R |^( io))j does not become zero in the number of intervals used to fit 
the staircase function to it. Inasmuch as h(o) is equal to the area 


under 



the value of h(o) 


can be made as small as desired by 


integrating farther out along the scale. The computations near t = 0 
usually require the most time and effort. 


Lightly Damped Airplane Response in Sideslip Angle 

In order to demonstrate the use of the incomplete time response to 
a unit Impulse on flight-test data, a typical pulse type of Input and 
response in sideslip of a large modern flexible swept -wing bomber is 
analyzed. Time histories of the side slip -angle response p(t) and 
rudder input 8 r (t) are shown in figure 5- The Incomplete response 0(t) 
(flight record cutoff at A = 5-4 seconds) and the complete time history 
of the Input S r (t) are used. In figure 6 , the Fourier transform ratio 
of the incomplete response 3(t) and Input S r (t) of the maneuver are 
shown. 


In order to compute the incomplete time response to a unit Impulse, 
it was first necessary to compute R J^-(ia>)J 

data of figure 6 . A staircase function was fitted to the curve of 


R 


=£-(M 

L 5 r J 


(shown in fig. 7 ) from the 

tted to the curve of 
of figure 7 by using a frequency Interval Ato = 0 .25 radian/sec. 


This staircase function was then used with table I of reference 2 to com- 
pute an Incomplete time response to a unit impulse shown in figure 8 . 


In order to check the accuracy of this incomplete time response to 
the unit impulse, the data of figure 8 were used with the original rudder 
pulse input 5 r (t), shown in figure 5 as the forcing function, to compute 

a time history of the sideslip angle p(t). The Duhamel method outlined 
In appendix A of reference 3 with a time interval At = 0.04 second was 
used for these calculations. A comparison of this computed 0(t) response 
with the measured 0 (t) response from the flight test is shown in the 
lower part of figure 5 by the circular symbols. This comparison Indicates • 
the magnitude of error involved in the total computation procedure 
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(transferring the data in time-history form to the response to a unit 
impulse and finally to the time response to an arbitrary input - by means 
of Duhamel' s integral . ) 

The Duhamel process was also carried out in the frequency plane. 

This computation was made by using tables I and III of reference 2. The 
process consisted of forming the Fourier transform of the arbitrary input 
and multiplying it by the ratio of the Fourier- transform of the incomplete 
response p(t) to the complete input & r (t), shown in figure 6, to give 

the Fourier transform of the desired response . The inverse of this 
Fourier transform then yields the desired time history of the response 
up to 90 percent of the cutoff time. The time response thus obtained Is 
shown in figure 5 by the square symbols. 


DISCUSSION 


In order to discuss the results presented, the manner of-estimating 
end corrections is first considered. The usual procedure is to know the 
mathematical relation describing the free oscillation of the system 
and to determine the parameters which fit into this relation from the 
time history of the response. For example, in figure 1 the relation 

describing the end correction is of the form Ce“ at sin(cut + 0) , where 
a is the damping parameter and C is a constant. These parameters 
can readily be determined by means of the logarithmic decrement, period, 
and phase of the time history of the response. Thus, for this simple 
case, an accurate frequency response can be obtained. It is important 
to note, however, In figure 2 that the contribution of the end correction 
to the frequency response can be considerable . Another Interesting pro- 
perty shown in figure 2 is that the cutoff time A appears in the fre- 


quency plane as the period 


2 rt 

A 


of an oscillation in both the amplitude- 


ratio and the phase-angle plots when an end correction is not applied. 

For example, in case 2 the period would be = 0.452 radian/sec. If 

13 • 9 


the application of an erroneous end correction permits such oscillations 
to persist, they could be misinterpreted as structural modes in the 
frequency -response analysis of data obtained from a flexible airplane . 


If-no end correction is added and the incomplete time response to 
the unit Impulse Is computed, the Incomplete time response i6 the same 
as the complete-system time response to the unit impulse up to 90 percent 
of the cutoff time • A. (See fig. 4.) This condition means that the 
incomplete, time response to the unit impulse is capable of completely 
defining the system from zero time up to A without the need of an 
end correction. 
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For many actual cases, such as the one shown in figure 5 > , the esti- 
mation of a valid end Correction is difficult, even though its mathemati- 
cal form may he derived from the equations of motion. The difficulty 
arises from attempting to determine the parameters which fit into this 
relation (for example, determining the damping from p(t) hy using the 
logarithmic-decrement technique for both the high- and the low-frequency 
inodes) . For such a system, the incomplete time response to the unit 
impulse shown in figure 8 has several advantages . It utilizes the data 
available without need of an. end correction and accurately defines the 
system in the time plane from zero time up to within 90 percent of the 
cutoff time. The incomplete time response to the unit impulse can be 
used in the Duhamel integral to compute the response to any arbitrary 
input up to the cutoff point. 

It has been assumed that the response to the unit impulse exists for 
the system in question; this assumption is usually the case. A necessary 
tho ugh not sufficient condition for the existence of the inverse Fourier 
transforms is that 


lim RjH(ioi)J = 0 

CD —foo 


In order to perform some of the computations shown in this paper it 
was necessary to extend the tables of reference 2. 


CONCLUDING REMARKS 


In the' present investigation it was found that the incomplete time 
response of a lightly damped dynamic system to a unit impulse could be 
accurately determined up to at least 90 percent of the cutoff time of 
the response from a knowledge of the time history of the complete input 
together with the incomplete time history of the response. In addition, 
this incomplete time response to a unit impulse may be used to predict, 
by means of Duhamel 1 s integral, the incomplete time response of the 
system to an arbitrary input up to at least 90 percent of the cutoff time. 
Alternatively, the Duhamel process may also be carried out in the fre- 
quency plane by using the arbitrary input and incomplete frequency response. 

The advantage of the incomplete time response to a unit impulse is 
that this method utilizes the data available and does not require an 
estimated end correction in order to describe adequately high-order 
lightly damped linear systems . Erroneous estimates of the end correction 
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made in obtaining the complete frequency response of a lightly damped 
linear system may cause a false impression of structural modes to appear 
in the frequency response of the system. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 17 , 1956. 
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Figure 2 .- Fourier transform ratio 5 of the system for three cutoff points. 
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Figure 2.- Concluded. 
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Figure 3 .- Real part of the Fourier transform ratio ^ of the system for 
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three cutoff points. 
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Figure 4.- Time response to a unit impulse of the system for three 

cutoff points. 
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Sideslip 



Time, t, sec 

Figure 5.- Typical Incomplete time history of a response in sideslip 
angle due to rudder deflection for a modem flexible swept -wing 
airplane. 
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Figure 8.- Incomplete time response to a unit impulse as obtained from 
the data of figure 5 (X - 5-4 sec). 
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